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CN i ABSTRACT 



In order to constrain the physical processes that regulate and downsize the Active 
Galactic Nucleus (AGN) population, the predictions of the MOdel for the Rise of 
[ GAlaxies aNd Active nuclei (morgana) are compared to luminosity functions (LFs) 

of AGNs in the optical, soft X-ray and hard X-ray bands, to the local BH-bulge 
mass relation, and to the observed X-ray number counts and background. We also 
give predictions on the accretion rate of AGNs in units of the Eddington rate and 
• on the BH-bulge relation expected at high redshift. We find that it is possible to 

\^ I reproduce the downsizing of AGNs within the hierarchical ACDM cosmogony and that 

. the most likely responsible for this downsizing is the stellar kinetic feedback that arises 

in star-forming bulges as a consequence of the high level of turbulence and leads to a 
massive removal of cold gas in small elliptical galaxies. At the same time, to obtain 
good fits to the number of bright quasars we need to require that quasar-triggered 
O ' galactic winds self-limit the accretion onto BHs; however, the very high degree of 

complexity of the physics of these winds, coupled with our poor understanding of 
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^ ' steepens considerably with respect to the observed one at bulge masses < 10^^ Mq; this 



■ it, hampers more robust conclusions. In all cases, the predicted BH-bulge relation 

> 



problem is related to a known excess in the predicted number of small bulges, common 
to most similar models, so that the reproduction of the correct number of faint AGNs 
. , is done at the cost of underestimating their BH masses. This highlights an insufficient 

^ . downsizing of elliptical galaxies, and hints for another feedback mechanism able to act 

on the compact discs that form and soon merge at high redshift. The results of this 
paper reinforce the need for direct investigations of the feedback mechanisms in active 
galaxies, that will be possible with the next generation of astronomical telescopes from 
sub-mm to X-rays. 

Key v^rords: galaxies: formation - galaxies: evolution - galaxies: ISM - quasars: general 



1 INTRODUCTION 

In the last years evidence has grown of the importance of the 
interaction between AGNs and their host galaxies. This in- 
teraction is now deemed to be of fundamental relevance for 
understanding the assembly not only of the supermassive 
Black Holes (BH) responsible for AGN activity, but also of 
the galaxies, in particular their spheroidal components (el- 
lipticals or spiral bulges). Observations of nearby galaxies 
show the existence of a well-defined correlation between the 
mass or central velocity dispersion of bulges and the mass 
of the hosted BHs (Kormendy & Richstone, 1995; Magor- 
rian et al., 1998; more recent determinations are found, e.g., 
in Marconi & Hunt, 2003 and Haring & Rix, 2004). More- 
over the observed mass function of these BHs is consistent 



with that inferred from quasar^ luminosities, under simple 
assumptions for the radiative efficiency and accretion rate 
(Soltan 1982; Cavaliere & Padovani 1988; Salucci et al. 1999; 
Yu & Tremaine 2002; Shankar et al. 2004; Marconi et al. 
2004; see also Haiman, Ciotti & Ostriker 2004). On the other 
hand, observations at high redshift highlight that quasars 
and radio-loud AGNs are hosted in elliptical galaxies (see 
e.g. Dunlop et al. 2003). 

An indirect evidence of the BH - bulge connection relies 
in the parallel evidence of downsizing of both populations. 
The evolution of the LFs of AGNs in the soft and hard 

^ In this paper we refer to luminous AGNs as quasars, with no 
reference to their radio emission. 
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X-ray bands reveals that the number density of fainter ob- 
jects peaks at a lower redshift with respect to brighter ones 
(see i.e. Ueda et al. 2003; Hasinger, Miyaji & Schmidt 2005; 
Barger et al. 2005; La Franca et al. 2005). This evidence 
of downsizing is often referred to as the "anti-hierarchical" 
behavior of BH growth, and is confirmed by the analysis of 
Merloni (2004) , based on the fundamental plane of accreting 
BHs (Merloni, Heinz & di Matteo, 2003), and by Marconi et 
al. (2004). Recently, the GOODS collaboration determined 
the LF of low- luminosity quasars at z ~ 4 (Cristiani et al. 
2004; Fontanot et al. 2006a), revealing a dearth of objects 
with respect to naive extrapolations of Pure Luminosity 
Evolution models anchored to the SDSS bright quasars (Fan 
et al. 2003). This evidence, confirmed by the COMBO 17 sur- 
vey at higher luminosities (Wolf et al. 2003) is again in line 
with the downsizing trend mentioned above. 

A similar trend has longly been claimed by some au- 
thors to be observed in galaxies (see, e.g., Cowie et al. 1996; 
Treu et al. 2005; Papovich et al. 2006; Fontana et al. 2006), 
especially in ellipticals, whose stellar mean ages, metal en- 
richment and level of alpha-enhancement are known to cor- 
relate positively with the stellar mass (see, e.g. Thomas 
1999; Renzini 2004). These evidences point to a formation 
scenario characterized by a quick burst of star formation 
rapidly followed by a strong wind which expels the resid- 
ual Inter-Stellar Medium (ISM); this wind should take place 
earlier in more massive galaxies (Matteucci 1994) and may 
be triggered by the quasar shining (Granato et al. 2001). In 
this framework it is also possible to explain the high metal- 
licity of quasar hosts (Matteucci & Padovani 1993; Hamann 
& Ferland 1999; Romano et al. 2002; D'Odorico et al. 2005). 

The influence of AGNs may not be limited to galaxies: 
jets from radio galaxies are now one of the most promis- 
ing candidates for quenching cooling flows in galaxy clusters 
(McNamara et al. 2005; Voit & Donahue 2005; Fabian et al. 
2006; see also the simulations of Quilis, Bower, & Balogh 
2001; Dalla Vecchia et al. 2004; Ruszkowski, Bruggen & 
Begelman 2004; Zanni et al. 2005; Sijacki & Springel 2006), 
and thus limit the mass of the most massive ellipticals (Ben- 
son et al. 2003; Bower et al. 2006; Croton et al. 2006). 

The problem of the demography of the AGN popula- 
tion in the context of galaxy formation in ACDM models 
has been addressed by many authors (see, e.g., Kauffmann 
& Haehnelt 2000; Monaco, Salucci & Danese 2000; Cav- 
aliere & Vittorini 2000, 2002; Cattaneo 2001; Granato et 
al. 2001, 2004; Mahmood, Devriendt & Silk 2004; Bromley, 
Somerville & Fabian 2004; Menci et al. 2003, 2004, 2006; 
Cattaneo et al. 2005; Vittorini et al. 2005; Bower et al. 2006; 
Croton et al. 2006; Hopkins et al. 2006; Malbon et al. 2006). 
In particular, the downsizing of galaxies and AGNs is not 
straightforward to reproduce in the hierarchical cosmogonic 
scenario, where larger dark matter (DM) halos form on av- 
erage later by the merging of less massive halos. As far as 
AGNs are concerned, Monaco, Salucci & Danese (2000) and 
later Granato et al. (2001; 2004) suggested that a delay in 
the shining of quasars in small bulges, motivated by stellar 
feedback, could explain this trend. Unfortunately, a more 
precise assessment of the success of the models listed above 
in predicting the downsizing of AGNs is not easy in most 
cases, as only a few of these papers (like Menci et al. 2004) 
explicitly compare to data in terms of the number density 
of AGNs as a function of luminosity. 



It is worth noting that in these papers different ideas 
have been proposed on the origin of the BH - bulge relation. 
Several proposed that the feedback from the AGN is able to 
self-limit the masses of both spheroids and BHs, (see, e.g., 
Ciotti & Ostriker 1997; Silk & Rees 1998; Haehneh, Natara- 
jan & Rees 1998; Fabian 1999; Granato et al. 2004; Murray, 
Quataert & Thompson 2005; Monaco & Fontanot 2005). 
Such self- limit at ion is supported also by the results of N- 
body simulations (Di Matteo et al. 2003, 2005; Kazantzidis 
et al. 2005; Springel, Di Matteo & Hernquist 2005). How- 
ever, the other authors successfully reproduced the BH - 
bulge relation simply assuming a proportionality between 
star formation (in bulges) and accretion rates, implicitly 
determined by the mechanism responsible for the almost 
complete loss of angular momentum of the gas accreting on 
the BH . The difference between the two approaches (self- 
limitation or proportionality between star formation and ac- 
cretion) was recently clarified by Monaco & Fontanot (2005). 

This paper is the second of a series devoted to describe 
the MORGANA code for the formation and evolution of galax- 
ies and AGNs. The code is described in full detail in Monaco, 
Fontanot & Taffoni (2006; hereafter paper I), while the abil- 
ity of the model to reproduce the early formation and as- 
sembly of massive galaxies is demonstrated in Fontana et 
al. (2006) and Fontanot et al. (2006b). Here we focus on 
the self-consistent modeling of accreting BHs and on their 
feedback on the host galaxies. By reproducing - or failing to 
reproduce - the observed properties of the AGN population, 
in particular the LFs in the optical, soft and hard X-ray 
bands, and the BH - bulge correlation at z = 0, we obtain 
valuable constraints on the physical processes involved. We 
also demonstrate our ability to roughly reproduce the ob- 
served X-ray counts in the soft (0.5-2 keV) and hard (2-10 
keV) bands and background from 0.2 to 300 keV, and give 
predictions on the BH - bulge relation at high redshift and 
on the Eddington ratios of the emitting objects. 

The paper is organized as follows. In section |5| we de- 
scribe the main properties of the MORGANA model, the ac- 
cretion onto BHs and its feedback on the forming galaxy, the 
runs used in the paper, and the procedure adopted to com- 
pute the AGN LFs, number counts and X-ray background. 
In section 13 we compare our results with available data and 
give some predictions on AGN and BH properties. The re- 
sults are discussed in section |1] while section |S] gives the 
conclusions. Throughout this work we assume the concor- 
dance ACDM cosmological model with Q.a = 0.7, Q.m = 0.3, 
Ho = 70 km s"^ Mpc-^ and cts = 0.9. 



2 MODELING AGNS IN morgana 

2.1 MORGANA 

MORGANA has been presented preliminary by Monaco & 
Fontanot (2005), and is described in full length in paper 
I. Here we give just a brief outline of the main physical pro- 
cesses included and discuss in some detail a few relevant 
points. 

The MORGANA model follows the typical scheme of semi- 
analytic models, with some important differences. Each DM 
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halo contains one galaxy for each progenitor^; the galaxy 
associated with the main progenitor is the central galaxy. 
Baryons in a DM halo are divided into three components, 
namely a halo, a bulge and a disc. Each component contains 
three phases, namely cold gas, hot gas and stars. For each 
component the code follows the evolution of its mass, metal 
content, thermal energy of the hot phase and kinetic energy 
of the cold phase. 

The main processes included in the model are the fol- 
lowing. 

(i) The merger trees of DM halos are obtained using 
the PINOCCHIO tool (Monaco et al. 2002; Monaco, Theuns 
& Taffoni 2002; Taffoni, Monaco & Theuns 2002). 

(ii) After a merging of DM halos, dynamical friction, 
tidal stripping and tidal shocks on the satellite (the smaller 
DM halo, with its galaxy at the core) lead to a merger with 
the central galaxy or to tidal destruction as described by 
Taffoni et al. (2003). 

(iii) The evolution of the baryonic components is per- 
formed by numerically integrating a system of equations for 
all the mass, energy and metal flows. 

(iv) The intergalactic medium infalling on a DM halo 
is shock-heated, as well as the hot halo component of merg- 
ing satellites (which is given to the main halo) and that of 
the main halo in case of major merger (Msat > 0.2 x A/tot). 
Following Wu, Fabian & Nulsen (2001), shock-heating is im- 
plemented by assigning to the infalling gas a specific ther- 
mal energy equal to 1.2 times the specific virial energy, 
—0.5Uii/Mh (where Mh and Uu are the mass and binding 
energy of the DM halo). 

(v) The profile of the hot halo gas is computed at each 
time-step by solving the equation for hydrostatic equilib- 
rium with a polytropic equation of state and an assumed 
polytropic index 7p = 1.2. 

(vi) Cooling of the hot halo phase is computed by in- 
tegrating the contribution of each spherical shell, taking 
into account radiative cooling and heating from (stellar and 
AGN) feedback from the central galaxy; if Tgo and pgo are 
the temperature and density of the hot halo gas extrap- 
olated to r = 0, fihotmp its mean molecular weight and 
Ts — rn/Cnfw the scale radius of the halo (of radius rn and 
concentration Cnfw), then the mass cooling fiow Mco,h re- 
sults: 



M,„,H = ^^xX(2/(7, 



^cool.O 



1)) 



where 



^cool.O 



3kTg{r cool) fihotmp 

2pgo(Acool — Fhcat) 



(1) 



(2) 



Here Acooi is the metal-dependent Sutherland & Dopita 
(1993) cooling function and the heating term Fhoat is com- 
puted assuming that the energy flow i5hw,H fed back from 



^ Each DM halo forms through the merging of many halos of 
smaller mass, called progenitors. At each merging the largest halo 
survives (it retains its identity), the others become substructure 
of the largest one. The main progenitor is the one that survives all 
the mergings. The mass resolution of the box used for computing 
the merger trees sets the smallest progenitor mass, as explained 
in section [2. 51 



the galaxy (including SNe and AGN) is given to the cooling 
shell: 



-Ehw,H 



47rriJ(2/(7p - 1)) 1 Pgo 



phot mp 



(3) 



Clearly these equations are valid if Fhcat < Acooi, otherwise 
no cooling flow is present. In equations and |3 the integral 
X(a) is defined as J^'"'"/^ {l-a[l-ln{l + t)/t]}°'t'^dt, with 

a = [3rvi,(7p - l)clZ{l + Cnfw)]/{7pTgo[(l + c^fw) ln(l + 
Cnfw) — Cnfw]} (Tvir being the virial temperature of the halo). 
The cooling radius rcooi is treated as a dynamical variable 
whose evolution takes into account the hot gas injected by 
the central galaxy (A/hw.n): 



A'/, 



'-'co,H 



47rp9(rcooi)r2^oi 



(4) 



(vii) The cooling gas flows into the cold halo gas phase; 
this is let infall on the central galaxy on a dynamical time- 
scale (computed at Tcooi)- This gas is divided between disc 
and bulge according to the fraction of the disc that lies 
within the half-mass radius of the bulge. In case of a disc- 
less bulge the disc size is estimated as Ru = O.TArn, where 
A is the spin parameter of the DM halo; a more precise com- 
putation is performed as explained in the next point when 
the disc accumulates a signiflcant amount of mass. 

(viii) The gas infalling on the disc keeps its angular mo- 
mentum; disc sizes are computed with an extension of the 
Mo, Mao & White (1998) model that includes the contribu- 
tion of the bulge to the disc rotation curve. 

(ix) Disc instabilities and major mergers of galaxies lead 
to the formation of bulges. We also take into account a pos- 
sible disc instability driven by feedback; this is explained in 
section IT^ In minor mergers the satellite mass is given to 
the bulge component of the larger galaxy. 

(x) Star formation and feedback in bulges and discs are 
inserted following the model of Monaco (2004a). For discs 
of gas surface density Ecom.d and fraction of cold gas /com.d 
the timescale for star formation t*.D is: 



t*,D = 9.1 



Scolc 



1 Mo pc- 



/ /cold,D Y 
V 0.1 ) 



Gyr 



(5) 



Due to the correlation of /coM.d and EcoM.d (galaxies with 
higher gas surface density consume more gas), this relation 
is compatible with the Schmidt law. For bulges the straight- 
forward Schmidt law is used: 



--coldjE 



1 M© pc- 



Gyr 



(6) 



In both cases, hot gas is ejected to the halo (in a hot galactic 
wind) at a rate equal to the star-formation rate (as predicted 
by Monaco 2004a), though massive bulges with circular ve- 
locity Vb > 300 km s"^ are able to bind the T ~ lO'' K hot 
phase component. The thermal energy of this re-heated gas 
is not scaled to the DM halo virial energy but is set equal 
to the energy of exploding SNe (assumed to be 10^^ erg) 
times an efficiency which is left as a free parameter. The 
best-fit value for the efficiency, 0.7, is very similar to the 
0.8 value suggested by Monaco (2QQ4b), who estimated the 
energy lost in the destruction of the host molecular cloud. 

(xi) In star-forming bulges cold gas is ejected in a cold 
galactic wind by kinetic feedback due to the predicted high 



4 Fontanot et al. 



level of turbulence driven by SNe. This is described in sec- 
tion 

(xii) When the hot halo phase is heated beyond the 
virial temperature, it can leave the DM halo in a galactic 
super- wind, at a rate: 



Mho 



(7) 



where iJhot.H and iSvir are respectively the actual and 
virial thermal energy of the hot phase (-Bvir /A/hot, h = 
— 0.5 (7h/A/h), Mhot.H its mass, taound its sound-crossing 
time and /wind is a free parameter set to the value of 2. A 
similar thing happens to the cold halo gas when it is acceler- 
ated by stellar feedback. To compute the time at which the 
ejected gas falls back into a DM halo (and is shock-heated) 
the merger history of the DM halo is scrolled forward in 
time until a halo is met with circular velocity larger than 
the (sound or kinetic) velocity of the gas at the ejection 
time. 

(xiii) Metal enrichment is self-consistently modeled in 
the instantaneous recycling approximation. 



2.2 Stellar feedback 

According to the model of Monaco (2004a), the regime of 
stellar feedback in a galaxy depends mainly on the density 
and vertical scale- length of the galactic system. In thin sys- 
tems, like spiral discs, the SNe exploding within a single 
molecular cloud give rise to super-bubbles that quickly blow 
out of the system, so that only a small fraction of their en- 
ergy is injected into the ISM, while most energy is given to 
the halo (~ 80 per cent according to Monaco 2004b). This is 
called adiabatic blow-out regime, and results in an ISM with 
thermal pressure of ~ 10^ K cm and a level of turbulence 
quantified by a velocity dispersion of clouds of ~ 7 km s~^ 
(see paper I). In thick systems, characterized by a large ver- 
tical scale-length or a high surface density, the super-bubbles 
do not manage to blow out of the system, which is then in 
the so-called adiabatic confinement regime, where most SN 
energy, both thermal and kinetic, is injected into the ISM. 
This results in higher thermal pressure and velocity disper- 
sion of clouds, well in excess of the ~ 6 km s~^ value found 
in spiral discs (see, e.g. Kennicutt 1989). 

A reference set of parameters for stellar feedback is in- 
troduced in paper I. Most parameters are fixed by requiring 
to reproduce the properties of local galaxies, or do influ- 
ence the predictions of the AGN population in a modest or 
rather predictable way, so we concentrate on varying only 
two mostly relevant parameters, leaving the others fixed to 
their standard values. The first of these parameters is re- 
lated to the amount of kinetic feedback in thick systems 
(i.e. in bulges; see point (x) of the list given in section ITTt . 
As shown in paper I in an equilibrium condition where the 
injection of kinetic energy from SNe is counter-balanced by 
the dissipation by turbulence, the velocity dispersion of cold 
clouds CTcoid scales with the star- formation timescale t* as: 



Ccold 



1 Gyr 



-1/3 



km s 



(8) 



The normalization parameter cro depends on many uncer- 
tain details, like the driving scale of turbulence. In thick 



systems, due to the high efficiency of (both thermal and ki- 
netic) energy injection, ctq is likely to be higher than in thin 
systems; this is observationally demonstrated by Dib, Bell 
& Burkert (2006). In conjunction with the much lower star- 
formation time-scale, this can lead to significant values of 
Ccoid in bulges. The resulting cold wind Mcw,b is quantified 
as: 



= A'/c,B-Pur 



funb 



(9) 



where Afc.B is the bulge cold gas mass, Punb the probability 
that a cold cloud with Maxwellian velocity distribution and 
rms velocity (Tcom overtakes the escape velocity of a bulge 
with circular velocity Vb and Vunb is the average velocity 
of the unbound clouds. We will show in the following that 
kinetic feedback in bulges plays a very important role in 
limiting faint AGNs at high redshift. We will test ao values 
ranging from to 90 km s"'^. 

The second parameter of stellar feedback that is consid- 
ered here is Eum, the threshold gas surface density for the 
switch to the thick system regime. This can be explained 
as follows: as a result of the strong cooling flows at high 
redshift, and of the assumption that the cooled gas settles 
on a disc, high-redshift discs may have very high surface 
densities of cold gas, sufficient to let them switch to the 
adiabatic confinement regime, typical of star-forming thick 
systems (bulges) and characterized by a higher velocity dis- 
persion of clouds. This process is beautifully seen in some 
high-redshift starburst galaxies with velocity fields typical of 
rotating discs and remarkably high gas velocity dispersions 
(of order of ~ 50 km s~^; Forster Schreiber et al. 2006). In 
these conditions the transport of angular momentum within 
the gaseous disc is more efficient, so that these objects are 
very likely to eventually evolve into bulges. This mechanism 
is implemented in a very simple way by stimulating a bar 
instability (which amounts to moving half of the disc mass 
to the bulge) whenever the gas surface density of the disc 
overtakes a value Eum; actually, radial fiows typical of bars 
are observed in some of the Forster Schreiber et al.'s (2006) 
objects, so this implementation could be realistic. The ref- 
erence value for the Sum parameter is suggested by Monaco 
(2004a; see also paper I) to be 300 M0 pc~^; the observa- 
tions mentioned above suggest that this may be a conser- 
vative choice, but varying the parameter in the range from 
100 to 500 Mq pc^'^ does not lead to very different results. 



2.3 Accretion onto BHs 

This part of the code is described in a simplified way in 
paper I and with slight differences in Monaco & Fontanot 
(2005). A seed BH of 10^ M0 is assigned to each DM halo, 
irrespective of its mass (see, e.g., Volonteri, Haardt & Madau 
2003 for a justification). Gas can accrete onto the BH only 
after having lost nearly all of its angular momentum J; fol- 
lowing Granato et al. (2004) we assume that this low-J gas 
accumulates in a reservoir of mass Mresv , from which it can 
accrete onto the BH. The first step in the loss of angular mo- 
mentum is connected to the same processes that lead to the 
formation of bulges; as a consequence, only the bulge cold 
phase can flow in the reservoir. Further losses of J may be 
connected to turbulence, magnetic fields or radiation drag 
(Umemura 2001); all these mechanisms are driven by star 
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formation, so the rate of accumulation of gas into the reser- 
voir, Miowj, will be related to the bulge star-formation rate 
Msi,B- In the simplest case the two flows will be propor- 
tional: Miowj = /lowj A^af.B"^- A more general relation be- 
tween Miowj and Msi,B is obtained assuming a power-law 
dependence with exponent Qiowj between the two quanti- 
ties: 



Ml, 



/lowjMsf.B 



100 Mq yr- 



(10) 



For aiowj = 1 this relation is equivalent to that of Granato 
et al. (2004), motivated by the radiation drag mechanism of 
Umemura (2001), while for Ofiowj 7^ 1 the /lowj parameter is 
scaled to a reference star-formation rate of 100 M© yr~^. For 
instance, the loss of angular momentum triggered by cloud 
encounters will likely have aiowj = 2. A similar approach is 
used by Cattaneo et al. (2005) 

The gas in the low-J reservoir accretes onto the BH at 
a rate determined by the viscosity of the accretion disc; this 
accretion rate Mvisc is found by Granato et al. (2004) and 
in paper I, but we report it here for completeness: 



ai f Aires. \ , Mbh \ 



(11) 



where (tb is the ID velocity dispersion of the bulge and fcaccr 
is theoretically estimated by the authors to be ~ 0.001. Ac- 
cretion is limited by the Eddington rate MeH/iEd, where 
t-Ed — 4 X 10^ yr is the Eddington-Salpeter timescale. The 
resulting system of equations for the BH mass Mbh and 
reservoir masse Mrcsv is (see also paper I): 



Mb 



Mr, 



= min Mv 



Mbh \ 



M, 



Mb 



(12) 



2.4 AGN feedback and quasar-triggered winds 

AGN activity releases a huge amount of energy so that, al- 
though the mechanisms for transferring it into the ISM are 
not very clear, this energy may easily trigger a massive galac- 
tic wind, able to remove all ISM from the galaxy. This would 
mark the end of the star formation episode that, according 
to the model described above, stimulated the accretion onto 
the BH. The details of the onset of such winds are very un- 
clear, so we decide to insert winds in the model in two ways. 

In both cases we use the same criterion for triggering 
the wind. This is motivated in Monaco & Fontanot (2005) as 
follows; the UV-X radiation of the AGN is able to evaporate 
some 50 Mq of cold gas for each Mq of accreted mass. When 
this evaporation rate overtakes the star-formation rate by a 
factor of order unity (0.3 in that paper), then the effect of 
the AGN radiation is sufficient to influence the ISM in a 
significant way. Removing one parameter (the evaporation 
efficiency and the triggering parameter are degenerate), the 
triggering condition can be written as: 



Mbh > /qwMst,E 



was named kr, 



(13) 



The quasar-wind parameter /qw determines the critical ac- 
cretion rate above which the wind is triggered, and takes 
values of order of 10"^ - 10"^; Monaco & Fontanot (2005), 
on the basis of the theoretically estimated 50 Mq of evap- 
orated gas per Mq of accreted matter, used paxameter val- 
ues corresponding to /qw = 0.006 to obtain reasonable self- 
regulated BH masses; we thus use 0.006 as a reasonable ref- 
erence value. Clearly, the criterion of equation 1131 is very 
similar to equation 1101 for aiowj = 1, which however refers 
to the build-up of the reservoir, not to the accretion rate of 
the BH. A modeling of the delay between loss of angular mo- 
mentum and accretion onto the BH, dictated by eguation llll 
is then necessary to use the criterion of equation 1131 

A massive removal of cold gas can take place only if 
the AGN is powerful enough to perform the work. Such a 
self-regulating mechanism is able by itself to produce a BH 
- bulge relation compatible with the one observed at z = 0; 
following again Monaco & Fontanot (2005) we put a second 
condition for the triggering of the wind, requiring that the 
mass of cold gas Mc,b to be removed from a bulge of mass 
A/b is not too large: 



M,. 



Mb 



< 0.21 



Mb 



4 Mq yr-i 



Mb 



10" Mq 



(14) 



A third condition for the triggering of a wind is set by 
requiring the BH to accrete in a radiatively efficient way; to 
achieve this we require that the accretion rate is more than 
1 per cent of the Eddington rate: 



Mbh > 0.01 



Mb 



tEd 



(15) 



in Monaco & Fontanot 2005. 



This is motivated by the low radiative efficiency of the flow 
in case of low accretion. These three trigger conditions are 
those proposed in the context of the Monaco & Fontanot 
(2005) model for the triggering of the wind, but their validity 
is wide enough to be used in a more general context. 

The modeling of the wind follows two routes. As a first 
option, "drying winds" are assumed to remove all the ISM 
from the bulge ejecting it to the halo. This is what happens 
if the wind is generated by an injection of kinetic energy 
coming directly from the accreting BH. The ejected gas is 
assumed to be heated to the inverse-Compton temperature 
of the AGN, T ~ 2 x lO'' K. Further accretion is possible 
from the reservoir, which is not depleted by the wind. This 
is motivated by the need to have a bright quasar phase after 
the wind has removed all the ISM, and may be justified by 
a bi-polar outflow that, after piercing the reservoir, gener- 
ates a blast wave that becomes symmetric while propagating 
in the ISM of the inner galaxy. As a second option, follow- 
ing the proposal of Monaco & Fontanot (2005), "accreting 
winds" are assumed to trigger further accretion onto the BH. 
This is what happens if the wind is generated throughout the 
galaxy by SNe^, then pushed away by radiation pressure of 
the shining AGN; in this case a part of the ISM is expected 
to be compressed to the center, so that a fraction /centre of 

* In the Monaco & Fontanot (2005) model tlie perturbation in- 
duced by runaway radiative heating of cold gas due to the shining 
of the quasar is able to trigger a change in the feedback regime 
that leads to the creation of a galaxy-wide outflowing cold shell, 
which is then effectively pushed out of the galaxy by radiation 
pressure. 
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the residual cold bulge gas is given to the reservoir of the 
BH. The reference value for this parameter is set to 0.002, 
the case in which ~ 20 per cent of the gas is compressed 
to the centre (as suggested for instance by Mori, Ferrara & 
Madau 2002 in a different context), and 1 per cent of this gas 
is able to lose its angular momentum and flow to the reser- 
voir; the results are rather insensitive to the precise value of 
this parameter. For simplicity we neglect any star formation 
connected to this compressed gas; this is done to minimize 
the effect of the highly uncertain mechanism of quasar winds 
on the host galaxy. Clearly, the case /centre = corresponds 
to the drying wind. 

Another feedback process that takes place, mostly when 
BHs are accreting in a radiatively inefhcient regime, is the 
heating of the hot halo gas by AGN jets, which can quench 
the cooling flows in large DM halos at low redshift. The 
importance of this process, which is complementary to the 
quasar winds discussed above, is highlighted also in Fontanot 
et al. (2006b). As explained in detail in paper I, we incor- 
porate a self-consistent implementation of this feedback by 
injecting the energy from the accreting BH to the hot halo 
gas each time the accretion rate is less than 1 per cent of Ed- 
dington; in this case the radiative efficiency in jets is known 
to be highest (see, e.g., Merloni, Heinz & di Matteo 2003). 
For higher accretion rates, we let only 10 per cent of the 
emitted energy heat the hot halo gas, on the ground that 
~10 per cent of quasars are radio-loud. A remarkable result 
highlighted in paper I (see its appendix B) is that the ef- 
fectiveness of the quenching depends on the mass resolution 
of the tree, so that particle masses < 10^ Mq are required. 
However, a proper sampling of the quasar population re- 
quires volumes well in excess of lO'^ comoving Mpc^; even 
for 512'^ particles this can be achieved only with particle 
masses larger than that limit. We then resort to a simpler 
recipe for the quenching of the cooling flow, with a philos- 
ophy similar to that of Bower et al. (2006) and Croton et 
al. (2006) (see also Cattaneo et al. 2005 and Kang et al. 
2006): we estimate the energy that the BH accreting at the 
highest radiatively inefficient rate (i.e. at 0.01 times the Ed- 
dington ratio) would give to the hot halo gas, then quench 
the cooling flow each time this energy is higher than that 
lost by cooling. As demonstrated in paper I, this "forced 
quenching" is able to mimic (with a slight overestimate) the 
effect of the self-consistent quenching described above even 
for poor mass resolutions. The mass acquired by the BH in 
quenching the cooling flow does not contribute to the AGN 
LFs, because the accretion is supposed to be radiatively inef- 
ficient by construction, and is not considered in the BH final 
mass. We stress that this "forced quenching" procedure is 
used in all the models presented below. 

2.5 Runs 

We use the same pinocchio runs introduced in paper I, in 
particular two 512^ pinocchio realizations of boxes of 200 
and 150 comoving Mpc {h — 0.7), with the cosmological 
parameters given in the Introduction. The mass particles 
are 2.4 x 10^ and 1.0 x 10® M©, so the smallest halo con- 
tains 50 particles, for a mass of 1.2 X 10" and 5.1 x 10^° 
Mq, while the branches of the DM halo merger trees start 
at a mass of 10 particles, corresponding to 2.4 x 10^" and 
1.0 x 10^° Mq. To sample rare objects like bright quasars a 



large box is required; however, as shown in paper I, there is 
some difference in the results when passing from the 200 to 
the 150 Mpc box, most notably in the amount of star forma- 
tion at z > 4 and in the effectiveness of the self-consistent 
quenching of the cooling flow (see above, section ETHl . which 
is poor for the larger box. We then decide to use the large 
box with the forced quenching procedure described above, 
which must be considered as a numerical trick to mimic the 
result obtained at higher mass resolution. As a matter of 
fact, the forced quenching procedure is slightly more effec- 
tive than the self-consistent one. In the following we will 
present only the results obtained with the 200 Mpc box, but 
we have checked that very similar results (with poorer sam- 
pling of the bright end of the LFs) are obtained with the 
higher-resolution box and physical quenching; the details of 
the quenching procedure have a very modest influence only 
on the high- luminosity tail of AGNs at z < 1. 

Considering then the 200 Mpc box, the stellar mass 
of the typical galaxy contained in the smallest DM halo at 
2: = is ~ 10® Mq; this is assumed as the completeness limit 
for the stellar mass function. For each run we compute the 
evolution of (up to) 100 trees (i.e. DM halos at z = 0) per 
logarithmic bin of halo mass of width 0.5 dex. This implies 
that while all the most massive halos are considered, smaller 
halos are randomly sparse-sampled. To properly reconstruct 
the statistical properties of galaxies we assign to each tree a 
weight Wtree equal to the inverse of the fraction of selected 
DM halos in the mass bin. 

The simulated comoving volume of 8 x 10^ Mpc'^ sets 
an upper limit to the number density of objects that 
can be studied with sufficient statistics. This limit depends 
on the probability of seeing an accretion event (with a given 
duty cycle) at a given redshift z in the box, and is computed 
as: 



where the limit refers to 10 objects in the box, tBd - the 
Eddington-Salpeter time - is used as a fiducial duration of 
an accretion event, and fbox is the cosmological time spanned 
by the box at the redshift z. For a box length of 200 Mpc 
this function takes values ranging from 10~^ Mpc~^ at z = 
to 2 X 10~* at z = 5. This way we cannot sample the 
brightest quasars, characterized by bolometric luminosities 
well in excess of lO*'^ erg s~^. This allows us to address the 
bulk of AGN activity at < z < 5, but not to consider the 
important problem of the assembly of bright quasars at very 
high redshift z ~ 6. This topic will be addressed elsewhere, 
using a larger simulated volume. 

2.6 Computing LFs and X-ray background 

As mentioned in paper I, the information on galaxies is out- 
put on a time grid of 0.1 Gyr. Due to the short duty cycle 
of the BH accretion events, this grid is too coarse to sam- 
ple properly the AGN activity. Then, information on BH 
accretion for all the galaxies is given at each integration 
time-step, whenever this accretion is significant; we use a 
limit Mbh > 1.76 x 10~^ Mq yr~^, corresponding to bolo- 
metric luminosities in excess of lO**^ erg s~^ . With this limit, 
the detailed information is issued only for a small fraction 
of integration time-steps. In particular, MORGANA outputs 
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the cosmological time t, the integration time-step At, the 
BH mass Mbh, the mass of the reservoir Mrosv, the bulge 
total mass Mb and cold gas mass Mc,b, the accretion rate 
onto the BH Mbh- Typically, each accretion event spans 
many contiguous timostops; wo treat each time-step as an 
independent event with duty cycle At. The accretion rate 
is converted into bolometric luminosity Lboio through the 
equation Lboio = ?7Mbh, where we assume that the accreted 
mass is converted into radiation with an efficiency of ?7 = 0.1. 
Each event is then counted Wtree times to correct for the 
sampling of the merger trees. 

We choose to test our model against observed LFs in the 
B, soft X-ray (0.5-2 kcV) and hard X-ray (2-10 kcV) bands. 
The hard X-ray band is the most useful one to compare with; 
its main advantage lies in the low level of extinction suffered 
by the radiation, especially for objects at high redshift for 
which an even harder radiation (in rest-frame terms) is ob- 
served. In this band we can assume that most objects are 
visible, with the only exception of Compton-thick AGNs, 
characterized by a hydrogen column density of Nh > lO^** 
cm~^; these however may be significant contributors of the 
X-ray background (see, e.g.. La FYanca et al. 2005). More- 
over, as the observed hardness ratio of a source allows one 
to estimate Nh, it is possible to estimate the unabsorbed 
hard-X flux. We compare our predictions with analytic fits 
of the observed LF, corrected for absorption, proposed by 
Ueda et al. (2003), Bargcr et al. (2005) and La Franca et al. 
(2005). 

Unfortunately, data at such high energies are rather 
sparse, so it is useful to compare also with the better sampled 
soft X-ray band. The best available data (Miyaji, Hasinger 
& Schmidt 2001; Hasinger, Miyaji & Schmidt 2005) are re- 
stricted to unabsorbed objects with Nh < W^^ cm~'^, so 
we need to correct for the fraction of absorbed objects be- 
fore comparing model and data. One possible way would be 
to assign an Nh to each AGN event, then selecting only 
those that satisfy the selection criterion; this however would 
decrease the statistics of the model LF. We prefer then to 
compute the fraction of unabsorbed objects in luminosity 
bins, given the Nh distribution, and correct our LF by that 
fraction. For the Nh distribution we use the luminosity- and 
redshift-dependent one proposed by La Franca et al. (2005) . 

We also compare our model to the B-band LFs (Ken- 
nefick, Djorgovski & Meylan 1996; Fan et al. 2003; Groom et 
al., 2004; Wolf et al. 2003; Cristiani et al. 2004; Fontanot et 
al. 2006a) , that arc measured with the best statistics and in 
the widest redshift range. As a matter of fact, the most strin- 
gent constraint comes from the high-redshift (3.5 < z < 5.2), 
low-luminosity (Mb ~ -22 to -24) AGNs observed by 
GOODS (Gristiani et al. 2004; Fontanot et al. 2006a). To 
compute our predicted B-band LF of type-I objects we as- 
sume that the type I fraction depends on luminosity accord- 
ing to the correlation found by Simpson (2005). 

To transform from bolometric to band luminosities we 
use the Elvis et al. (1994) bolometric correction for the B- 
band, assuming a value of 10.4±2.0 for the ratio i^44oo/i',44oo 
and the bolometric luminosity Lboio- In the X-ray bands we 
adopt the bolometric corrections proposed by Marconi et al. 
(2004): 

Log{^^^) = 1.64 + 0.22Li2 + O.OI2L12 - 0.0015L?2 



Log{^^) = 1.53 -I- 0.24Li2 -|- 0.012L?2 - 0.0015Lf2 (17) 

where L12 = I/0(;(Lboio) — 12. Marconi et al. (2004) also pro- 
pose a luminosity-dependent bolometric correction for the 
B-band which is in agreement with the Elvis et al. correc- 
tion in the I/boio range of our interest. 

As a consistency check, we compare our models also to 
number counts and the X-ray background. To compute our 
predictions for these quantities we build a library of tem- 
plate spectral energy distributions (SEDs) as follows (see 
also Monaco & Fontanot 2005; Ballo et al. 2006). In the op- 
tical we consider the quasar template spectrum of Gristiani 
& Vio (1990) down to 538 A, extrapolated to 300 A using 
oc v-^ ''^ (following Risaliti & Elvis 2005). At shorter 
wavelengths (between 0.01 and 30 A) we use a power-law 
SED with a photon index of F = —1.8 and an exponen- 
tial cutoff exp(—_E/200keV) (see i.e. La Franca et al., 2005). 
The relative normalization between the optical-UV and the 
X-ray branches of the spectrum is constrained through the 
quantity aox (Zamorani et al., 1981): 

OLox = — -, T- (lOj 

lO3(l'2500/i^lkev) 

We use for aox the bolometric luminosity-dependent of aox 
value proposed by Vignali, Brandt & Schneider (2003): 

otox = -0.11 logiUisoo) + 1-85 (19) 

The interpolation between 30 and 300 A follows Kriss et al. 
(1999). We then produce a library of template spectra in the 
range Log(Z/boio) = [42.0, 47.5] (divided into bins of 0.1 dex 
in luminosity). For each accretion event in our morgana 
output we associate a template spectrum and an Nh value, 
extracted from the La Franca et al. (2005) distribution, 
which includes also Compton-thick objects (for which the X- 
ray flux is set to zero). The Nh absorption is computed using 
the Morrison & McCammon (1983) cross section. We also 
compute absorption by the ISM following Madau, Haardt & 
Rccs (1999); this is important only at the lowest energies. 
The integration in redshift is easily performed with the MOR- 
GANA output, as this spans the whole range of cosmological 
times from recombination to the present. This is at variance 
with Fontanot et al. (2006b), where the computation of the 
galaxy SEDs is much more demanding in terms of computer 
time, so that a very careful sampling of model galaxies is 
needed. 



2.7 Parameter space and models 

Starting from the standard set of parameters defined in 
paper I, we have investigated a limited subset of param- 
eters that influence significantly the AGN population but 
are varied in a limited range so as to give very similar re- 
sults in terms of galaxies. These are the galaxy feedback 
parameters ctq (which sets the level of kinetic feedback) and 
Slim (which sets the feedback-induced bar instability), and 
the quasar wind parameters /qw (which sets the trigger for 
quasar winds), /lowj (which regulates the quantity of gas 
flowing into the reservoir), aiow.T (which sots the scaling of 
the angular momentum loss with the bulge star-formation 
rate) and /centre (which sets the fraction of bulge gas that is 
available for accretion after the wind is triggered). 
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Model (TO Slim /qw /low.1 CtiowJ /centre 

STD 60 km s"! oo Mq pc-^ 0.003 1 

DW 60 km s-i oo Mq pc'^ 0.006 0.03 2 

AW 60 km s-i 300 Mq pc-^ 0.006 0.01 1 0.002 



Table 1. Parameter values for the models used in this paper. 

We have investigated this parameter space, finding more 
that one acceptable solution. Here we present the results ob- 
tained with three models, namely the standard model pre- 
sented in paper I (called STD), a model with drying winds 
and aiowj = 2 (DW) and a model with accreting winds 
Slim ~ 300 M© pc"'^ (AW). The relative parameters are 
reported in tabled It is worth noticing that the two mod- 
els with quasar winds have higher /iow.i values: in the STD 
model /lowj sets the the level of the BH - bulge relation, 
while in the wind models DW and AW the final BH mass 
is not set by /lowj but by the self-regulating action of the 
wind; in this case much higher /lowj values can be allowed, 
and this parameter regulates the early bulding of the mas- 
sive BHs and then the bright end of the AGN LFs. 



3 RESULTS 

In this section we compare the results of our model to obser- 
vations of LFs and number counts in the hard-X, soft-X and 
optical bands, and to the statistics of remnant BHs at z = 0. 
To best illustrate the constraints that can be obtained on the 
physical processes described above we show results for the 
three combination of parameters given in tableland intro- 
duced in section im The STD model refers to the standard 
choice of parameters presented in paper I, with the addition 
of the forced quenching procedure described in section VTM 
In this case quasar winds are not active (/q„ = 0). The re- 
sulting LFs for this model are always too steep (figures^^J 
though the z = Q BH - bulge relation is roughly reproduced 
(figure for Mb > lO" Mq. Within the range of param- 
eters allowed by the constraints of galactic observables (see 
paper I and Fontanot et al. 2006b) we have found no way 
to have shallower LFs. On the other hand, the introduction 
of quasar winds allows us to improve the agreement with 
AGN data without influencing much the results on galaxies. 
The DW model includes drying winds as follows. In order 
to flatten the AGN LF the Ofiowj exponent is set to 2, so 
that accretion is increased in the strong starbursts that give 
rise to big spheroids, and depressed in small bulges; in order 
to avoid a dramatic steepening of the BH - bulge relation 
a,t z = Q, drying winds are introduced to limit the mass of 
the most massive BHs. In this case we use /iow,t = 0.03 in 
place of the 0.003 value of STD. Another good combination 
of parameters (model AW) is found by allowing accreting 
winds (with /qw — 0.006, /lowj = 0.01 and /centre = 0.002) 
and setting Eum = 300 Mq pc~^ (model AW); in this case 
both the starbursts induced in discs with a high gas surface 
density and the secondary accreting episodes stimulated by 
the quasar wind increase the activity of bright quasars even 
with aiowJ = 1, 

Figures [HOI and El show that the DW and AW models 
flt nicely (though maybe not in detail) the observed LFs of 
quasars in all the bands considered. 



Figure 31 shows the same results in a different way. The 
predicted number density of AGNs in bins of bolometric lu- 
minosity is compared to the range of values inferred from 
observations (using mainly the analytic flt of the hard-X 
LFs of Ueda et al. 2003, Barger et al. 2003 and La Franca et 
al. 2005 for z < 3.5, and the results of Fontanot et al. 2006a 
at z > 3.5). The three models are all shown. It is clear that 
the DW and AW models reproduce nicely the evolution of 
the AGN population, with the number distributions peak- 
ing at lower redshift for lower luminosity. In other words, the 
downsizing or anti-hierarchical behavior of AGNs can be re- 
covered in the context of the hierarchical ACDM cosmogony 
once feedback is properly modeled. 

Similar results are obtained with the 150 Mpc box, 
though the statistical limit for the number density of ob- 
jects (equation ^Jl is higher by a factor of 2.4, the decrease 
of the AGN activity at low redshift is marginally slower and 
the slight excess of faint AGNs at z ~ 2 is a bit larger. 
Clearly, some of these discrepancies can be absorbed by a 
further tuning of the parameters, which we deem worthless 
in this analysis. 

From a much more detailed analysis of the parameter 
space we infer that the physical process at the origin of the 
downsizing of AGNs is the kinetic feedback active in star- 
forming bulges. This is shown in flgure |51 where the DW 
model is shown with ao values ranging from to 90 km 
s~^ (very similar results are obtained with the other mod- 
els). Kinetic feedback causes a strong ejection of cold gas 
in small bulges, thus decreasing the number of faint AGNs 
without changing much the number of bright quasars. From 
this comparison we obtain a best-flt value for ao of 60 km 
s~^; this is the most effective way to constrain this param- 
eter. It is also worth noticing in figure |S] the fundamental 
importance of the constraint coming from the abundance of 
high-z AGNs in the GOODS fields (Cristiani et al. 2004; 
Fontanot et al. 2006a). 

Figure |S] shows the BH - bulge relation at 2: = pre- 
dicted by the models, while figure Q shows the resulting 
mass function of BHs. All models, included STD, roughly 
reproduce the observed range of values for massive bulges 
{Mb > 10" Mq), with some modest overestimate by the 
wind models (especially AW), while a steepening is predicted 
at smaller masses. This figure allows to make some impor- 
tant points. First, despite its fundamental importance in 
demonstrating the connection between AGNs and their host 
bulges, this relation is not after all a very strong constraint, 
as the three models give similar fits (for Mb > 10^^ Mq) 
despite their remarkably different performance in terms of 
LFs and number densities of objects. Second, the scatter of 
this relation is remarkably low in the STD model, and this 
suggests that the recipes given in section ITU are too simple 
and some source of scatter is needed. Intriguingly, quasar 
winds give roughly the correct amount of scatter, though 
they tend to give too many massive BHs at z = 0. Third, 
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Figure 1. 2-10 keV X-ray LF. The lines are tlie predictions of the STD model (black solid line), DW model (red dashed line), AW model 
(blue dot-dashed line). The shaded region represents the Ueda et al (2003) estimate. The cyan horizontal bands give the limit below 
which the LF is affected by poor statistics in our run (equation I16i 



BHs in small bulges are too small compared to observations. 
This very important point will be deepened below. 



ground synthesis that uses similar ingredients as their pa- 
per. This agreement shows that also the population of faint 
sources is roughly reproduced by the wind models. 



3.1 X-ray number counts and background 

In order to strengthen the constraints on the AGN popu- 
lation we compare our model predictions to X-ray number 
counts (figure ISJ in the soft and hard X-ray bands. Apart 
from some possible overestimate at the brightest fluxes, the 
fit to the data of the two wind models (DW and AW) is 
very good, while the STD model overestimates the counts 
by roughly a factor of three, and this is due to the overesti- 
mate of the faint end of the LFs. 

Figure|5|shows the prediction for the X-ray background 
from 0.5 to 300 keV. The background predicted by the two 
wind models follows nicely the observed one, though the 
peak at ~ 30 keV is underpredicted by the DW model. 
This is in line with the results of La Franca et al. (2005), 
who claim a missing (though not dominant) population of 
Compton-thick AGNs, which is not included in this back- 



3.2 More predictions 

In order to quantify the amount of evolution in the BH-bulge 
relation we compare our model predictions with the results 
of Peng et al. (2006), who measured the evolution of the BH 
- bulge relation by defining the quantity r(2): 



Mb Mb 



(20) 



Mbh Mbh 

In figure irni we present the average evolution of F pre- 
dicted by the three models, but only for the galaxies with 
Mb > 2 X 10^^, which safely lie on the relation. The model 
STD show some degree of evolution, which is due to the fact 
that, consistently with Croton (2006), the fraction of stars 
formed in bulges is higher at higher redshift, so that a higher 
fraction of gas can accrete on the BH. On the other hand. 
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Figure 2. 0.5-2 keV X-ray LF. Lines refer to models as in figure^ In all panels the empty circles and squares refer respectively to the 
observations of Miyaji et al. (2001) and Hasinger et al. (2005). 



the wind models show a more pronounced increase of F with 
redshift, in better agreement with Peng et al. (2006). This 
was anticipated in Monaco & Fontanot (2005), who showed 
that if more mass is allowed to flow onto the BH until the 
winds limits further accretion, then the mass assembly of 
massive BHs is anticipated. This can help in reconciling the 
presence of bright quasars at high redshift, when very mas- 
sive galaxies were not yet assembled. 

Regarding the mechanisms responsible for the assem- 
bly of BHs, we find that the bulk of BH mass is gained from 
accretion. Figure fTTI shows the fraction of mass acquired by 
accretion (the rest is acquired by mergers) as a function of 
BH mass for the BHs found in z = galaxies, according 
to the three models. We notice that merging is significant 
for the most massive BHs, which reside in the most massive 
galaxies that experience complex merger histories (see, e.g., 
De Lucia et al. 2006; De Lucia & Blaizot 2006). We also 
notice that quasar winds, especially the accreting winds, in- 
crease the amount of mass acquired by accretion. 

Another important prediction is the average accretion 
rate of BHs in units of the Eddington rate. This quantity is 



needed to relate the accretion history of AGNs, estimated 
from the LFs, to the local BH mass function. As each accre- 
tion event is subdivided by our model into many sub-events, 
one per integration interval, we show this quantity averaged 
over bins of bolometric luminosity and redshift and weighted 
by the width of the time bin. Figure fT?l shows the results for 
the STD, DW and AW models. While the STD rates are 
very low, especially for the highest accretion rates, the Ed- 
dington ratios of the other models are rather high at high 
redshift and decrease at low redshift, especially for low AGN 
luminosities. This prediction is compared with the observa- 
tional results of KoUmeier et al. (2006) , based on a sample 
of 407 AGNs from the AEGIS survey, with BH masses esti- 
mated through a combination of line widths and continuum 
luminosities; the agreement is reasonably good, though Ed- 
dington ratios tend to be higher than the data at high red- 
shift. This shows that the reason why the STD model fails to 
produce bright quasars relies in the low accretion rates stim- 
ulated by star formation, while the main effect of winds is to 
force massive BHs to accrete at high rates. A similar predic- 
tion was given in Kauffmann & Haehnelt (2000), based on a 
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Figure 3. B-band LF. Lines refer to models as in figure [Tl Data are taken from Kennefick et al. (1996), Wolf et al. (2003), Croom et al. 
(2004), Fontanot et al. (2006a), as indicated in the figure. 



parametric description of the accretion timescale. Our more 
refined modeling allows us to propose that this behaviour 
can be obtained only by inserting quasar-triggered galaxy 
winds. Finally, the predicted dependence of the Eddington 
ratios on redshift and bolometric luminosity can be useful 
to relate the accretion history of BHs to the observed mass 
function of remnant BHs at 2 = 0. 



4 DISCUSSION 

The main message of this paper is that it is possible to in- 
clude accreting BHs in a complete galaxy formation model 
like MORGANA so as to reproduce the main constraints on the 
AGN population from the optical to the hard X-ray. How- 
ever, this good agreement is reached at the cost of meddling 
with the very uncertain physics of quasar-triggered galac- 
tic winds. By presenting two different good solutions, both 
based on motivated and plausible choices on how to insert 
quasar winds, we intend to stress that the parameter space 
is so wide, even in the restricted version presented here, that 



it is not possible to constrain in a unique way the complex 
processes at play simply by reproducing the statistical prop- 
erties of the AGN population. 

The second message of the paper is that, notwithstand- 
ing their non-uniqueness, the only acceptable solutions we 
find are based on quasar winds. As clear in figure 1121 the 
role of quasar winds is to force a high accretion rate for the 
massive BHs and reproduce the trend of increasing Edding- 
ton ratio with bolometric luminosity, in better agreement 
with data. Given the uncertainties on the physics of these 
events, it is possible to use this result as a suggestive clue 
on the role of winds but not as a strong evidence of their 
necessity. This caution is confirmed by the fact that most 
papers cited in the Introduction are able to reproduce bright 
quasars without advocating any such mechanism. However, 
there is a deep difference between our model, where the ac- 
cumulation of low- J gas is connected to star formation, and 
most other models where some cold gas is produced at a 
triggering event (mergers if not flybies) and accreted at an 
arbitrarily specified fraction of the Eddington rate. In this 
regard our model can be directly compared only to that of 
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Figure 4. Evolution of the QSO number density with redshift in different bins of bolometric luminosity, estimated from X-ray (Ueda et 
al. 2003; Barger et al. 2003; La Franca et al. 2005) and optical (SDSS, Fan et al. 2003; GOODS, Fontanot et al. 2006a) surveys. Lines 
refer to models as in figure [Tl 



Granato et al. (2004), who proposed a very similar model- 
ing of the accretion onto the BH, but used a lower value of 
their parameter equivalent to /lowj; in their case the wind 
takes place when the star-formation episode has consumed 
nearly all the gas, and their winds do not limit strongly the 
BH mass. Their model should then be more similar to our 
unsuccessful STD model. But the real difference between 
their approach and MORGANA relies in the treatment of star 
formation and feedback; in particular, they assume that the 
formation of an elliptical takes place in a single episode and 
that the cooling gas goes to the bulge, while in our model 
and in the same context the cooling gas would settle on a 
disc^ until it loses angular momentum through bar instabil- 
ities, mergers or feedback (see the discussion on the role of 
Slim in section [2.211 . As a consequence, in morgana many 



of the stars that end up in bulges are formed in discs and 
do not contribute to the loss of angular momentum of gas 
(this is the main reason for the scatter in the BH - bulge 
relation). Our higher values of /lowj allow to have stronger 
accretion events and higher BH masses at high redshift, but 
then a limiting mechanism, like quasar winds, is necessary 
to avoid to overshoot heavily the BH - bulge relation. In this 
regard, the modest overestimate of large BH masses by the 
DW and AW^ models could be fixed by a more careful study 
of the parameter space, but we find hints that winds are a 
key ingredient to achieve agreement with many observables. 

Our model is also comparable to that of Cattaneo et al. 
(2005), who relate the accretion rate to the star- formation 
rate in a similar way as equation 1101 with their equivalent 
of the aiowj parameter equal to 1.5. Unfortunately, their 



^ In this context our choice of letting the cooled gas flowr in the 
bulge is not of particular relevance, as it is important only when 
most of the galactic baryonic mass is in a bulge, so after the 
galaxy has formed. 



^ In the AW case, where the discrepancy is more pronounced, 
we have verified that this overestimate is due to many minor 
accretion events that take place at very low Eddington rates and 
thus do not contribute to the AGN activity. 
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Figure 5. Evolution of tlie QSO number density for the DW model as a function of ctq. The black solid line refers to ctq = km s~^, 
the magenta dotted line to (tq = 30 km s~^, the blue dashed line to ctq = 60 km s~^ (the best-fit value) , the red dot-dashed lines to 
iTo = 90 km s~^. 



relatively small box (100 Mpc//i with 256^ particles) does 
not allow to reach strong conclusion on the prediction of 
their model for bright quasars. 

The third message is that kinetic feedback in star- 
forming bulges is a very good candidate for the main mech- 
anism responsible for the downsizing of AGNs. This point 
however needs some discussion. An overall good agreement 
between models (DW or AW) and data is obtained only for 
Mb > 10" M0 and Mbh > 10* M©; at smaller values the 
BH - bulge relation steepens considerably, especially in the 
DW model. We think that this discrepancy points to an in- 
trinsic problem of MORGANA that is shared with most other 
galaxy formation models. As shown in paper I, the predicted 
stellar mass function of bulge-dominated galaxies does not 
have a broad peak at ~ 10^" Mq, as suggested by observa- 
tions, but presents a power-law tail of small objects which is 
just below that of discs; this problem is shared, for instance, 
by the model of Croton et al. (2006). Recently, Fontana et al. 
(2006), based on the GOODS-MUSIC sample (Grazian et al. 
2006), reconstructed the stellar mass function of galaxies up 



to z ~ 4 and compared their results with N-body and semi- 
analytic models including MORGANA. All the models were 
found to overpredict the number density of small galaxies at 
z ^ 1; many of these galaxies are the likely progenitors of 
small bulges. Similar results are obtained with the Garching 
model (De Lucia, private communication). The mechanism 
of kinetic feedback can suppress effectively star formation in 
small bulges, but, as mentioned above, cannot hamper stars 
that form in discs to get into bulges. As a consequence, the 
formation of BHs is suppressed but the number of bulges 
is not, and the BH - bulge relation steepens (remarkably, 
the steepening is slightly less evident in the AW case, where 
the discs with high gas surface density are transformed into 
bulges). In other words, a good fit of the number of AGNs 
and an overestimate of the number of bulges can only be 
compatible with a steep BH - bulge relation. 

The excess of small bulges is slightly over-corrected by 
kinetic feedback, in that the resulting BH mass function 
is low at small masses (figure Ell. This over-correction can 
be explained as follows: kinetic feedback, which is applied 
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Figure 6. Predicted Mb - J\/bh relation at 2 = for STD, DW 
and AW models. The thick diagonal line gives the observed mean 
relation according to Marconi & Hunt (2003), the dotted lines 
highlight the observed scatter of 0.3 dex. 
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Figure 7. BH mass function at z = for the models, (lines refer 
to models as in figure 0, compared to the region allowed by the 
Shankar et al. (2004) and Marconi et al. (2004) estimates (shaded 
area). 



Figure 8. Cumulative Source Number Counts in X-ray bands. 
Left Panel: soft (0.5-2 keV). Estimates from Hasinger et al., (1993) 
{dotted region); Brandt et al., (2001) (empty star); Bauer et al., 
(2004) (empty circle); Hasinger (1998) (filled triangle); Zamorani, 
Mignoli & Hasinger (1999) (filled circle); Rosati et al. (2002) 
(empty squares). Right Panel: hard (2-10 keV). Estimates from 
Giommi et al., (2001) (empty triangles); Ogasaka et al. (1998) 
(filled square); Brandt et al., (2001) (empty star); Bauer et al., 
(2004) (empty circle); Rosati et al. (2002) (green empty squares). 
Lines refer to models as in figure [Tl 



at a rather high level, limits the accretion on small BHs 
(Mbh =i 10® - 10^ Mq), but these BHs are anyway hosted in 
larger bulges (figures 1^1 and lTITl , so they are likely to have ac- 
cess to larger amounts of gas; then, a given (small mass) BH 
can give rise to more (low-luminosity) AGN activity than 
desired. The constraint on the total accretion leads then to 
a further underestimate of BH masses. But if the small BHs 
are to provide the correct amount of AGN activity, they 
should accrete on average at a higher rate than observed. 
To test this idea we compare our predictions with the obser- 
vations of Ballo et al. (2006), who estimated stellar masses, 
bolometric luminosities and BH masses (based on the BH 
- bulge relation ai z = 0) for a hard X-ray-selected sample 
of faint AGNs hosted in galaxies at 0.4 < z < 1. Figure [T^ 
shows the Eddington fraction of the DW and AW models 
compared to the Ballo et al. (2006) data for Mbh > 1.6 x lO'^ 
M0 and for all BH masses. In particular, the thick lines show 
model predictions with the same mass limit; these are in 
very good agreement with the corresponding observational 
points (the larger triangles). This means that the accretion 
rate of these BHs is correctly reproduced by the model. This 
result does not depend on the details of the quenching of 
cooling fiows by AGN feedback, as the involved BH masses 
are rather small; indeed, identical results are obtained in 
the higher resolution box with the standard quenching pro- 
cedure. Besides, the accretion rate of all BHs (thin lines), 
which is dominated by the smaller objects, is tendentially 
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Figure 9. Reproduction of tlie X-ray background. Observations 
are highlighted by the shaded areas and are taken from Revnivt- 
sev et al. (2003), Worsley et al. (2004), Georgantopoulos et al 
(1996), Vecchi et al. (1999), Lumb et al., (2002), De Luca & 
Molendi (2004), Kinzer et al., (1997), Gruber et al. (1992). Lines 
refer to models as in figure [Tl 



Figure 11. Scatter-plot of the predicted fraction of BH mass 
acquired by accretion as a function of BH mass, for BHs in z = 
galaxies. The three panels give the results of the STD (upper), 
DW (mid) and AW (lower) model. 
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Figure 10. Evolution in redshift of the BH— bulge relation 
through the F parameter feQuation l20i for the three models. Lines 
refer to models as in figure Data points are taken from Peng 
et al. (2006). 



Figure 12. Predicted average accretion rates in units of Edding- 
ton; thick lines refer to models as in figure^ Thin lines represent 
the variance of the DW model, which is representative also for the 
other models. Stars refer to data from KoUmeier et al. (2006). 
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Figure 13. Predicted Eddington ratios as a function of tlie bolo- 
metric luminosity for BHs more massive than 10^-2 M0 (thick 
lines, red dashed for the DW model, blue dot-dashed for the AW 
model) and for all BHs (thin lines, styles as above). Data are from 
Ballo et al. (2006), large points refer to BHs more massive than 
lO''^'^ Mq, smaller points to smaller BHs. The thick lines should 
be compared with the large points, the thin lines with all the 
points. 



larger than observed (the thin lines must be compared with 
both large and small triangles). This trend is not very strong, 
and could be due to a lack of small BHs in the dataset; how- 
ever, as discussed in the Ballo et al. (2006) paper, the hard 
X-ray selection makes this possible bias weak if not absent. 
This comparison then confirms that, to compensate for the 
excess of small bulges, our model generates too few small 
BHs [M < 10^ M0) that, being hosted in relatively larger 
bulges, accrete at a higher rate. The kinetic feedback mech- 
anism remains then a very good candidate for downsizing 
AGNs, but the value proposed here of 60 km s~^ is proba- 
bly overestimated. 

The discrepancy described above points to some missing 
feedback mechanism active at high redshift able to downsize 
the small star-forming discs in the same way as kinetic feed- 
back does for the small bulges. Bulge formation by feedback, 
explained in section helps but does not solve the prob- 
lem, so other mechanisms are required. 

Other authors have obtained good agreement between 
the prediction of similar semi-analytic models and data in 
terms of the downsizing of the AGN population, so it is 
worth wondering if other mechanisms than kinetic feedback 
can give similarly good results. In the case of Granato et 
al. (2001, 2004), the downsizing is obtained by delaying the 
shining of quasars, as originally suggested by Monaco et al. 
(2000). This delay is justified by stellar feedback, but the 
model does not make a distinction between thermal and ki- 
netic feedback. Menci et al. (2003, 2004, 2006; see also Vit- 
torini et al. 2006) obtain similarly good results, although the 



level of predicted downsizing may be slightly less than that 
required by data (Menci, private communication). In their 
case the downsizing is obtained by stellar feedback (param- 
eterized as usual with a /3 coefficient which is a power-law of 
the disc velocity) and by the insertion of galaxy-galaxy in- 
teractions as a further trigger of AGN activity, a mechanism 
which penalizes the small satellites. Similarly, Cattaneo et 
al. (2005) with a standard feedback recipe and an accretion 
rate similar to eauation llOl (but without reservoir) obtain a 
roughly good match to the optical LF of quasars. In all these 
cases stellar feedback is treated in a simple and parametric 
way. As shown and discussed in Monaco (2004a) and Paper 
I, feedback in discs can lead to the reheating of an amount 
of cold gas roughly equal to the star-formation rate, which 
means /3 ~ 1, so the only way to have more cold gas ejected 
to the halo is to accelerate it more than heating it. We are 
then rough agreement with the papers mentioned above in 
stating that the suppression of faint AGNs at high redshift 
is due to stellar feedback, but our more refined modeling 
allows us to draw stronger conclusions on the details of the 
physical mechanism at play. 



5 CONCLUSIONS 

This paper belongs to a series devoted to presenting MOR- 
GANA, a new model for the joint formation and evolution of 
galaxies and AGNs, and is focused on the role of feedback in 
shaping the observed properties of AGNs and their relation 
with galaxies. Our analysis confirms that models based on 
the ACDM cosmogony are able to roughly reproduce at the 
same time both the properties of AGNs (presented here) and 
the properties of galaxy populations (presented in paper I, 
Fontana et al. 2006 and Fontanot et al. 2006b). 

This model has been extended, through a careful mod- 
eling of AGN SEDs, to predict many QSO properties, espe- 
cially in the X-rays. It reproduces nicely LFs and number 
counts in the soft and hard X-ray bands, and the measured 
background from 0.5 to 300 keV, with a possible underesti- 
mate of the peak which confirms the claim for a significant 
though not dominant population of Compton-thick sources. 
This implies that the bulk of AGN accretion is roughly re- 
produced. 

The agreement between model and data allows us to 
draw these conclusions. 

Quasar triggered winds are necessary in our model to 
reproduce the number density of bright quasars. The param- 
eter space of quasar winds is discouragingiy wide; within a 
very limited sub-space, characterized by "only" six param- 
eters, we find two good solutions, which is a clear sign that 
our knowledge on the physics of accretion onto BHs and 
their interaction with galaxies is still too poor to draw firm 
conclusions. In any case, the idea that quasar winds are a 
necessary ingredient of galaxy formation is worth pursuing. 
The two wind solutions that we propose are based on the 
same triggering criterion for the quasar wind, requiring that 

(i) the accretion rate is high enough to perturb the ISM, 

(ii) the gas mass to remove is not too high, (iii) the AGN is 
accreting in a radiatively efficient mode. The first solution is 
based on a tilt of the relation between bulge star-formation 
rate and loss rate of angular momentum, compensated by 
"drying winds", where the kinetic energy injected by the 
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central engine causes a complete removal of ISM from the 
host bulge. The second solution is based on the "accreting 
wind" mechanism proposed by Monaco & Fontanot (2005), 
where the wind is generated throughout the ISM, so that 
a fraction of the cold gas is compressed to the center and 
stimulates further accretion. A good fit of the AGN LFs is 
then obtained by assuming, as explained in paper I, that 
discs with high gas surface density lose angular momentum 
(and then become bulges) because of the expected change 
in the feedback regime (Monaco 2004a) and the consequent 
increase of the velocity dispersion of clouds (as observed by 
Forster Schreiber et al. 2006). 

We propose that kinetic feedback in star-forming bulges 
is a very good candidate for downsizing the AGN popula- 
tion. Indeed, the high velocity dispersion of cold gas in star- 
forming bulges can lead to a massive removal of ISM, lead- 
ing to a suppression of less luminous AGNs at high redshift, 
while larger elliptical galaxies can retain their gas more eas- 
ily. At lower redshift larger BHs may accrete at lower rates, 
giving rise to the low-luminosity population dominating the 
X-ray background. 

We make the following predictions: (i) the BH-bulge re- 
lation is already in place at high redshift, though high BH 
masses in relatively small bulges may be found at 2 > 2 if 
quasar winds are active; (ii) the Eddington ratios of accret- 
ing BHs depend both on redshift and bolometric luminosity, 
with low values for the faint AGNs at < z < 1 that give the 
bulk of the hard X-ray background (consistent with Ballo et 
al. 2006); (iii) most BH mass is acquired through accretion 
more than through mergers, especially if quasar-triggered 
winds are switched on. 

There is an important point of discrepancy between our 
model and data, in that the predicted BH - bulge relation at 
Mb < W^^ Mq, or Mbh < 10* Mq, is significantly steeper 
than observed; in the same mass range the predicted mass 
function of BHs is too low. To our understanding, this is 
connected to a known excess of small bulges with respect to 
observations (see paper I) and to the excess of small-mass 
galaxies predicted by most (N-body or semi-analytic) mod- 
els at z ~ 1 (Fontana et al. 2006). We have tested this idea 
by comparing our model predictions of the Eddington ratios 
of faint AGNs at 0.4 < z < 1 with the data of Ballo et al. 
(2006); while the agreement is excellent for M > 1.6 x 10^ 
M0, smaller BHs accrete at a higher rate than observed, 
and this compensates for the lack of smaller BHs. This dis- 
crepancy highlights that while kinetic feedback is efficient in 
downsizing AGNs by quenching the star formation in bulges, 
stars in small ellipticals may be born in discs, so that an- 
other mechanism is needed to produce the required level of 
downsizing of bulges. 

In conclusion, in developing this work we have gained 
interesting insight into the complex problem of the cosmo- 
logical rise of AGNs, highlighting two promising astrophysi- 
cal mechanism (kinetic feedback in star-forming bulges and 
quasar-triggered galaxy winds) and the need for a third, un- 
known one to improve the downsizing of elliptical galaxies. 
However, our poor understanding of the underlying physics 
hampers more robust conclusions. This demonstrates once 
again that the field of joint galaxy and AGN formation is 
observationally-driven and we need strategies to single out 
the physical mechanisms at play. In particular, the main de- 
grees of freedom in the theory are related to the connection 



between BH accretion and star formation, which can be con- 
strained by estimating accretion rates, star-formation rates 
and BH masses in active galaxies; to the quasar-triggered 
winds, which can be constrained by observing warm and 
cold absorbers or Lyman-a blobs associated to quasars; and 
finally to the nature of stellar feedback, which can be con- 
strained by detailed observations of starburst galaxies. The 
next generation of telescopes will provide suitable tools to 
assess these topics. 
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